/*
 *  main.cpp
 *  T3nsors2
 *
 *  Created by Michael Barriault on 10-07-05.
 *  Copyright 2010 University of Guelph. All rights reserved.
 *
 */


#include "MW.h"
#include <cstdlib>
#include <string>
#include <iostream>
using namespace t3;
using std::string;
using std::cout;
using std::endl;

int main(int argc, char* argv[]) {
	int n = 8;
	double K = 0.5;
	double epsilon = 0.05;
	double s = 0.25;
	double T = 1.0;
	for ( int i=1; i<argc; i++ ) switch ( argv[i][0] ) {
		case 'n':
			n = stoi(((string)argv[i]).erase(0,1));
			break;
		case 'e':
			epsilon = stof(((string)argv[i]).erase(0,1));
			break;
		case 'K':
			K = stof(((string)argv[i]).erase(0,1));
			break;
		case 's':
			s = stof(((string)argv[i]).erase(0,1));
			break;
		case 'T':
			T = stof(((string)argv[i]).erase(0,1));
			break;
	}
	domain* O = domain::init(3,n, -0.2,0.2,	-0.2,0.2,	-2.0,2.0);
	MW S(O, K, epsilon);
	S.Solve(T,s);
	O->release();
	return 0;
}

